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Abstract. We study cosmological models that contain sterile neutrinos with eV-range masses 
as suggested by reactor and short-baseline oscillation data. We confront these models with 
both precision cosmological data (probing the CMB decoupling epoch) and light-element 
abundances (probing the BBN epoch). In the minimal ACDM model, such sterile neutrinos 
are strongly disfavoured by current data because they contribute too much hot dark matter. 
However, if the cosmological framework is extended to include also additional relativistic 
degrees of freedom beyond the three standard neutrinos and the putative sterile neutrinos, 
then the hot dark matter constraint on the sterile states is considerably relaxed. A further 
improvement is achieved by allowing a dark energy equation of state parameter w < —1. 
While BBN strongly disfavours extra radiation beyond the assumed eV-mass sterile neutrino, 
this constraint can be circumvented by a small v e degeneracy. Any model containing eV-mass 
sterile neutrinos implies also strong modifications of other cosmological parameters. Notably, 
the inferred cold dark matter density can shift up by 20-75% relative to the standard ACDM 
value. 
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1 Introduction 

There is mounting evidence from reactor and short-baseline neutrino oscillation experiments 
suggesting the existence of one or two sterile neutrinos with mass splittings relative to the 
active flavours in the neighbourhood of Am 2 ~ 1 eV 2 and fairly large mixing parameters 
(see, e.g., [1] for a review and global interpretation, and [2-4] for other recent analyses). In 
the early universe, flavour oscillations would bring these sterile states into thermal equilib- 
rium prior to neutrino decoupling at T ~ 1 MeV, thereby increasing the relativistic energy 
density (see, e.g., [5] for a review). This increase, commonly parameterised in terms of 
an additional contribution to the relativistic neutrino degrees of freedom AT e g, significantly 
modifies big-bang nucleosynthesis (BBN) of light elements, the cosmic microwave background 
(CMB) anisotropies, and the formation of large-scale structures (LSS). By the same token, 
these eV-mass sterile neutrinos would later play the role of a non-negligible hot-dark matter 
component. 

The cosmological verdict on eV-mass sterile neutrinos is somewhat mixed. Current 
CMB and LSS observations show a consistent preference for additional relativistic degrees of 
freedom beyond the standard model expectation of N e ^ = 3.046, with low to moderate sta- 
tistical significance [6-11]. On the other hand, if this putative radiation excess is interpreted 
in terms of sterile neutrinos, the usual hot dark matter limits constrain the sterile masses to 
the sub-eV regime. For example, in a 3+1 scenario consisting of three essentially massless 
active flavours and one fully thermalised sterile flavour, the sterile mass is constrained to 
m s < 0.5 eV [7]. Thus, cosmology appears to disfavour the existence of one or two sterile 
states that have mass and mixing parameters favoured by a global analysis of the laboratory 
data. 

However, both the tentative evidence for extra radiation and these hot-dark matter 
bounds are based on the simplest ACDM framework; extended models may provide different 
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answers. For example, if the dark energy equation of state parameter w and curvature £1^ 
are used as fit parameters, the presence of eV-mass sterile neutrinos favours a w value much 
smaller than —1 [12] . Here, we take the view that it may be more natural to look for extensions 
to ACDM within the neutrino sector alone, for example in terms of additional radiation. 
(However, we note in passing that in certain dynamical dark energy models such as the mass 
varying neutrino scenario, a strongly modified neutrino sector is responsible for driving the 
dark energy evolution. In this sense, changing w may not be an entirely unnatural solution 
to the massive sterile neutrino problem.) Multiple sterile right-handed states may exist and 
some of them may even be lighter that the eV range. Moreover, the existence of sterile 
neutrinos also provides a means for the generation of large neutrino chemical potentials [13] , 
which may be exploited to circumvent standard BBN bounds on the number of relativistic 
degrees of freedom. 

The premise of our work is that sterile neutrino states with eV-range masses are real 
and that they are fully thermalised prior to neutrino decoupling. In addition, we allow for 
radiation beyond that due to the three standard flavours (assumed to be effectively massless) 
and beyond the eV-mass sterile neutrinos. We consider both data from precision cosmology 
and the constraints on iV e jj imposed by BBN. 

More specifically, in section 2 we investigate the impact of eV-mass sterile neutrinos on 
the inference of cosmological parameters from CMB and LSS observations. In section 3 we 
evaluate the constraints on additional relativistic species during the BBN epoch, using mea- 
surements of the primordial abundances of several light elements. In section 4 we summarise 
and discuss our results. 

2 CMB and LSS 
2.1 Cosmological data 

The CMB anisotropies and the LSS distribution are sensitive to both additional relativistic 
degrees of freedom and the mass scale of light, free-streaming particles. To explore the im- 
pact of eV-mass sterile neutrinos on precision cosmology, we use CMB anisotropy data and 
their accompanying likelihood routines from the WMAP 7- year data release [14], as well as 
the ACBAR [15], BICEP [16], and QuAD [17] experiments. In addition, we use the halo 
power spectrum extracted from the SDSS-DR7 luminous red galaxy sample [18], and type 
la supernova (SN) data from the Union-2 compilation [19]. Finally, we impose a constraint 
on the Hubble parameter based on the Hubble Space Telescope observations [20] . We refrain 
from using observational data pertaining to very small length scales, such as estimates of 
small-scale density fluctuation amplitudes from the Lyman-a forest and cluster abundance. 
While these measurements are in principle an extremely powerful tool for constraining neu- 
trino masses, they are also currently dominated by systematic uncertainties. We have tested 
explicitly that adding data from the Atacama Cosmology Telescope (ACT) [21] does not alter 
our results; The main role of ACT is to break parameter degeneracies in CMB-only analyses. 
These degeneracies are however already broken by the addition of LSS data. 

We construct allowed regions in cosmological parameter space using standard Bayesian 
inference techniques. All posterior probability density functions are sampled using the 
Markov Chain Monte Carlo (MCMC) code CosmoMC [22]. 
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Table 1. Definitions of the model parameters and the prior ranges imposed on them in our CMB+LSS 
likelihood analysis. Uniform priors are used for all listed parameters. 



Parameter 


Symbol 


Prior 


Baryon density 




0.005 -> 0.1 


Cold dark matter density 


Wcdm 


0.01 -> 0.99 


Hubble parameter 


h 


0.4 -> 1.0 


Amplitude of scalar spectrum @ k = 0.05 Mpc -1 


log[10 10 yl a ] 


2.7 4 


Scalar spectral index 


n s 


0.5 ->■ 1.5 


Optical depth to reionisation 


T 


0.01 -» 0.8 


Number of extra massless neutrino degrees of freedom 
Dark energy equation of state parameter 


AiV ml 

w 


0^5 
-2.5 -> -0.5 



2.2 Cosmological models 

We consider three basic cosmological frameworks in which we embed light sterile neutrinos 
and analyse their consequences. 

1. The ACDM class of models is defined by a flat spatial geometry and six free parameters, 
{^b ; Wcdm, h, A s , n s , t}. See table 1 for their definitions and prior ranges. Within this 
framework, we consider four possibilities in the neutrino sector: (a) 3 massless neutrinos 
(standard ACDM), (b) 3 massless+1 sterile (0 eV), (c) 3 massless+1 sterile (1 eV), and 
(d) 3 massless+1 sterile (2 eV). 

2. In the second class of models, which we dub ACDM+AA r , we include additional rela- 
tivists degrees of freedom beyond the 3+1 standard and sterile neutrinos. We consider 
three scenarios: (a) 3+AN m \ massless+1 sterile (0 eV), (b) 3 + AN m \ massless+1 sterile 
(1 eV), and (c) 3 + AN m \ massless+1 sterile (2 eV). Such models are especially sensitive 
to BBN constraints to be discussed in section 3. 

3. The wCDM+AN framework is a variant of ACDM+AiV, in which we extend the dark 
energy sector to include the possibility that w ^ — 1. 

2.3 Goodness-of-fit 

We consider first the goodness-of-fit, as quantified by the best-fit effective x 2 > of the cos- 
mological models described in the previous section. Here the effective x 2 is defined as 
Xgff = — 21n£ max , where £ max is the maximum likelihood of the data given the model. 
Table 2 summarises the best-fit xts f° r ^ese models and, where appropriate, the preferred 
values of AA^i and w, using the data sets described in section 2.1. Some comments are in 
order. 

1. Within the ACDM framework, a scenario with 3 massless neutrinos plus one fully 
thermalised massless sterile species offers a slightly better fit than standard ACDM 
(Ax 2 ff = —3.16). However, once the sterile states have mass, the quality of the fit 
deteriorates. The 1 eV scenario is already marginally worse than standard ACDM by 
Axl s = 4.20. The 2 eV model may be deemed unacceptable (Ax 2 ff = 21.41). 

2. The situation improves when we allow for additional radiation (ACDM+AA r ). For 
example, if the sterile neutrino has a mass of 1 eV, the best-fit XeS ^ s comparable to 
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that of standard ACDM, albeit at the cost of admitting A.N m \ ~ 1.5 additional massless 
degrees of freedom. For a 2 eV sterile mass, we find AN m \ ~ 2.6 and XeS = 12.8. This 
last result can be compared with the u>CDM+f2fc model of reference [12], for which 
A%^ ff = 12 assuming a lighter 1.33 eV sterile neutrino. Therefore, introducing extra 
radiation appears to be somewhat superior to modifying the dark energy sector at 
resolving the sterile mass conundrum. 

3. Even more improvement is available if, in addition, we allow the dark energy equation 
of state parameter w to differ from —1 (wCDM+AiV). In this class of models, we see 
that a scenario with one species of 1 eV sterile neutrinos in fact provides a better fit 
to the data than does standard ACDM, with &Xeg = — 0.78, at the expense of two 
additional free parameters. 



Table 2. Best-fit A^g relative to the standard ACDM framework for the models described in 
section 2.2, using the data sets of section 2.1. We also show the best-fit values and 95%-credible 
upper and lower limits on w c( | m , and, where appropriate, on AA m i and w. 



Framework 


Neutrino sector 




Axcff 


AAT ml 


w 




k>cdm 


ACDM 


3 massless 













n i qo+0.0036 




3 massless + 1 sterile (0 eV) 




-3.16 








0-1299t~ 




3 massless + 1 sterile (1 eV) 




4.20 








o.i398tS:88£ 




3 massless + 1 sterile (2 eV) 




21.41 








147S+0-0075 
U ' 14M -0.0064 


ACDM+AiV 


3+AiV m i massless + 1 sterile (0 


eV) 


-3.54 


o oi +i - 12 

u - ui -0.01 






n , OO+0.023 
U.loo_o 005 




3+AiV m i massless + 1 sterile (1 


eV) 


2.26 


1 40+ 111 






0.166ir f 7 




3+AiV m i massless + 1 sterile (2 


eV) 


12.82 


9 r 7 +1.24 
-0.59 






0.192«;?^ 


wCDM+AiV 


3+AiV m i massless + 1 sterile (0 


eV) 


-5.38 


u - uy -0.09 


r- uu -o.i2 


n , OO+0.032 
u - ±oz -0.006 




3+AiV m i massless + 1 sterile (1 


eV) 


-0.78 


-1 oq+161 
L - ZO -0.75 


-l.llj 


-0.18 
-0.21 


n 1R4+ 035 

U.1D4_ 015 




3+AiV m i massless + 1 sterile (2 


eV) 


7.80 


2.48tJ; 7 g 


— 1.17j 


-0.23 
-0.22 


1 QC+0 032 

U.198_ 019 



2.4 Effects on other cosmological parameters 

We have seen that precision cosmological observations can reasonably accommodate one fully 
thermalised species of massive sterile neutrinos if we allow also for additional massless degrees 
of freedom and/or a non-standard dark energy equation of state. An interesting consequence 
is that the preferred values of other free parameters of the model also shift accordingly. 

The most notable example in this regard is the cold dark matter density u C( i m . Figure 1 
illustrates the shift in w c d m as a function of the sterile neutrino mass within the ACDM+AA^ 
framework. Figure 2 is similar, but for the wCDM+AN models. See also table 2 for the 
best-fit values and credible regions. Clearly, the larger the sterile neutrino mass, the larger 
the preferred value of cj c d m - In the case of a 2 eV sterile neutrino, the upward shift in w c dm 
can be as large as 75% in the u^CDM+AA^ model, relative to the standard ACDM inferred 
value. This shift in the cold dark matter density can have importance consequences for, e.g., 
the SUSY dark matter parameter space. 

Another affected parameter is the scalar spectral index n s , whose preferred region widens 
in the presence of additional light species, as was also seen in our previous analysis [7]. 
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Figure 1. 2D marginal 68%- and 95%-credible regions in the (AN m \, w c d m )-plane for three 
ACDM+A7V class models containing one thermalised sterile species of mass m s = eV (solid), 
1 eV (dashed), and 2 eV (dotted). 

-0.6 
-0.8 



-1 
-1.2 



-1.4 



-1.6 

0.1 0.12 0.14 0.16 0.18 0.2 0.22 0.24 

^cdm 

Figure 2. 2D marginal 68%- and 95%-credible regions in the (u>, w c d m )-plane for three wCDM+AiV 
class models containing one thermalised sterile species of mass m s = eV (solid), 1 eV (dashed), and 
2 eV (dotted). 




3 Big bang nucleosynthesis 

While not sensitive to the masses of neutrinos, BBN has long been used to probe the radiation 
content of the Universe at temperatures of order 1 MeV [23-25]. In this section, we explore 
the implications of the latest primordial element abundance measurements on the sterile 
neutrino scenario. We consider a general BBN model with three free parameters: the baryon 
density w b > an extra N s effective sterile neutrino species on top of the usual three fully 
thermalised standard neutrinos, and, eventually we also allow for the presence of a neutrino 
chemical potential £. Table 3 summarises the prior ranges for these parameters. 

3.1 Analysis 

We infer constraints on the BBN parameters with a modified version of the Markov Chain 
Monte Carlo (MCMC) sampler CosmoMC [22]. 1 The PArthENoPE [26] code is employed to pre- 
calculate the primordial element abundances on a grid of points in (wb, N s , £)-space. For each 

1 The numerical module for implementing the BBN likelihood in CosmoMC is available from the corresponding 
author (Jan Hamann) upon request. 
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Table 3. Definitions of the model parameters and the prior ranges imposed on them in our BBN 
likelihood analysis. Uniform priors are used for all listed parameters. 



Parameter Symbol 


Prior 


Baryon density Uh 
Number of sterile neutrinos N s 
Neutrino chemical potential £ 


0.01 -> 0.035 

0^5 
-0.2 -> 0.3 




Figure 3. Primordial abundances of D, 4 He and 7 Li as functions of the cosmological parameters. 
Top left: dependence on lo^ for iV s = £ = 0. Top right: dependence on N s for Wb = 0.0225 and £ = 0. 
Bottom: dependence on £ for Wb = 0.0225 and N s = 0. 

set of parameter values the abundances are then obtained from the grid by 3-dimensional 
cubic spline interpolation. We illustrate the dependence of the light element abundances on 
the cosmological parameters on selected slices through parameter space in figure 3. 

3.2 Uncertainties of input data 

The output of PArthENoPE is subject to theoretical uncertainties from two principal sources: 
the nuclear reaction rates and the free neutron lifetime r n . The former induce a ±1.6% 
error for the D prediction, negligible error for 4 He, and ±8% for 7 Li [27]. We fold these 
uncertainties into our definitions of the corresponding likelihood functions in section 3.3. 

The main source of uncertainty for the predicted 4 He abundance is the free neutron 
lifetime. The Particle Data Group recommends a value of Tjf DG = 885.7 ± 0.8 s [28]. This, 
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however, appears to be in strong tension with recent measurements performed by Sere- 
brov et al. [29], who found = 878.5 ± 0.8 s, and Pichlmair el al. [30], who measured 
= 880.7 ± 1.8 s. Although a re-analysis of neutron lifetime data shows that previous es- 
timates may have been biased by ~ +6 s [31], we will nonetheless consider the two extreme 
values and r n PDG in the following, in order to illustrate the possible impact of this uncer- 
tainty on our inference. Going from to t pdg shifts the Y p prediction by ~ 0.6%, whereas 
[ 2 D/H] and [ 7 Li/H] change by ~ 0.3% which is negligible compared with the effects of 
other uncertainties (nuclear rates and astrophysical measurements). 

3.3 Primordial element abundance data and likelihood functions 
Deuterium 

The primordial deuterium abundance [D/H] can be inferred from measurements of the 
absorption of quasar light in high-redshift low-metallicity hydrogen clouds. Pettini et al. [32] 
derive a value of log [D/H] = —4.55 ± 0.03 from the observation of the spectra of seven 
quasars. Taking into account the theoretical uncertainty of 1.6%, the corresponding likelihood 
function is 

(log [D/H] +4.55) 2 
" 21 ^ = A 0.034* <") 

Helium-4 

The mass fraction Y of 4 He is determined by measuring helium and hydrogen emission lines 
in low-metallicity HII regions of dwarf galaxies. Its primordial value Y p can then be estimated 
by treating Y as a function of metallicity Z and extrapolating to Z = 0. A linear regression 
of data from seven high-quality objects with a careful treatment of systematic uncertainties 
yields Y p = 0.2573^Q QQgg, with the additional restriction to theoretically meaningful positive 
slopes dY/dZ > [33]. 2 We therefore model the likelihood function as 



(y p -0.2573) z 
0.0033 2 



, if Y p > 0.2573, 



0.0088 2 ' p 



if y D < 0.2573. 



Lithium-7 



Lithium abundances have been determined from the spectra of metal-poor dwarf stars. At 
low metallicities, the lithium abundance appears to be independent of metallicity [37]. How- 
ever, the measurements are subject to significant systematic uncertainties about the stars' 
temperatures. Additionally, the measured plateau value may not be representative of the 
primordial abundance: 7 Li can be generated by cosmic rays, or depleted through Population 
III stars or diffusion processes within the dwarf stars [38]. 

The Particle Data Group estimates an average of [ 7 Li/H] p = (1.7 ±0.06 ±0.44) x lO" 10 , 
which is significantly lower than the standard BBN expectation for realistic baryon densities 
- the well-known lithium problem [39]. Approximating the 8% theoretical uncertainty with 

2 It has been argued that a linear regression may not necessarily be realistic in view of certain models of 
early 4 He production by Pop III stars [34], as for instance suggested in references [35, 36]. However, even 
in these cases, the difference between the primordial and the observed values of Y will not exceed ~ 0.01, a 
possibility that is well covered by the lower limit on Y p in equation (3.2). 



-7- 



an absolute error of c™? 1 ~ 0.4 x 10 10 and adding the errors in quadrature, we arrive at the 
following likelihood function for lithium: 

(lQi°rLi/H] p -1.7) 2 
"21n£ 7Li = — 2 , (3.3) 

which will be considered only where appropriate. 
CMB+LSS prior on the baryon density 

In addition to primordial abundance measurements, we also consider a prior on the baryon 
density uy, from CMB+LSS data as an independent constraint. Since the bounds on are 
in principle model-dependent, using results from a fit to the standard ACDM model as a 
prior for extended models would be technically incorrect. Instead we use the baryon density 
inferred within the ACDM+Ag model, oj h = 0.02255 ± 0.00049, to define our prior 

91 f ^ K ~ 0.02255) 2 
- 2 In vrcMB+LSS K) oc q OQQ492 " ( 3 - 4 ) 

In practice, however, present CMB+LSS data are sufficiently sensitive to break any potential 
degeneracy between N s and u^, so that the constraints on show no appreciable variation 
with respect to the inclusion or otherwise of the parameter N s . 

3.4 BBN with sterile neutrinos 

We first consider a BBN scenario with N s effective sterile neutrino species. The additional 
radiation energy density contributed by the sterile neutrinos increases the expansion rate 
at BBN, leading to a higher freeze-out temperature Tf. As a consequence, the equilibrium 
neutron-to-proton ratio n/p = exp(— Am/T), where Am is the neutron-proton mass differ- 
ence, is larger at neutron freeze-out. The resulting Y p is larger because almost all neutrons 
end up in 4 He. Similarly, a larger neutron lifetime also increases the expected Y p . 

We show the constraints on the (u^, N s ) parameter space from our analysis in figure 4. 
The top two panels illustrate the well-known fact that the deuterium abundance constrains 
mainly the baryon density, whereas constraints on N s are mostly driven by the helium data. 
However, once we impose the CMB+LSS prior on o^,, even the deuterium data alone signif- 
icantly constrain N s . This constraint is of interest in view of the helium data already being 
systematics-limited. 

For a combined fit of the deuterium and the helium data, the best-fit value of N s is 
0.86, with a 95%-credible upper limit of N s < 1.26 (or 1.24 if the CMB+LSS prior on cjb is 
included). Using the larger neutron lifetime , we obtain slightly lower values: a best-fit 
of 0.73 and a 95% upper limit of N s < 1.14. The one-dimensional marginalised posterior 
probability densities for N s are shown in figure 5. As already emphasised in reference [34], 
there is no strong indication for N s > from BBN alone - N s = lies well within the 90%- 
credible interval in all cases - owing to the relatively weak lower limit on Y p . While one fully 
thermalised sterile neutrino species is slightly favoured over N s = 0, two fully thermalised 
sterile neutrino species are clearly incompatible with the data in this scenario. 

3.5 Degenerate BBN with sterile neutrinos 

If we imagine for a moment that there were two sterile neutrino species today, how could this 
be reconciled with the results of the previous section? One possibility would be incomplete 
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0.035 




Figure 4. 2D marginal 90%- and 99%-credible regions in the (wb, N s ) -plane, assuming the 
BBN+sterile scenario. Black lines denote results from elemental abundance measurements alone 
(solid lines assume r„ = 878.5 s, dotted lines 885.7 s). The red lines include a CMB+LSS prior on 
Wb- Top left: D data. Top right: 4 He data. Bottom: D+ 4 He data. 




Figure 5. ID posterior pdf for N s , marginalised over Wb, for the BBN+sterile scenario and D+ 4 He 
data. Same colour/line coding as in figure 4. 



thermalisation, such that the effective N s is smaller than 2; this scenario could be confirmed 
if future CMB+LSS data should find iV e g < 5. Alternatively, the sterile neutrinos could be 
the decay products of a heavy particle species between BBN and decoupling. In this case, 
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N s could be smaller than 2 at BBN, but equal to or larger than 2 at decoupling. 

A third possibility is the degenerate BBN scenario [40], in which all standard neutrinos 
share a common non-zero chemical potential £ [41-43]. The presence of a chemical potential 
affects BBN in two ways. Firstly, it contributes an additional term to the effective radiation 
density, 



45 r 
7 



2 (£/tt) 2 + (C/tt) 



(3.5) 



Secondly, the degeneracy in the electron neutrinos modifies the equilibrium neutron-to-proton 
ratio, 

( Am \ / s 

n/p = exp — — £ J . (3-6) 



For £ ~ O(0. 1), it is the latter effect that is most relevant for the resulting element abun- 
dances: a positive £ will reduce the number of available neutrons with respect to the £ = 
case and thus suppress Y p . This suppression can be used to circumvent the upper bound on 
N s found in the previous section. 




Figure 6. 2D marginal 90%- and 99%-credible regions in the (£, N s ) -plane marginalised over Wb, 
assuming the degenerate BBN+sterile scenario and r n = 878.5 s. The black lines denote results from 
D+ 4 He data, while the red lines also include a CMB+LSS prior on uj^,- 

Figure 6 shows the preferred region in the (£, iV s )-plane in the degenerate BBN scenario. 
In the 3-dimensional parameter space (u\,,N a , £), the deuterium and the helium data can only 
constrain two directions, leaving an unconstrained one which admits arbitrarily high values 
of N s . Only by adding the CMB+LSS prior on are we able to obtain any constraint at 
all: we find a best-fit of N s = 0.79 and a 95%-credible upper limit of N B < 2.56. Conversely, 
if we assume N s = 2, a positive non-zero chemical potential is required: 0.03 < £ < 0.14 
(95%-credible interval), with a best-fit value of £ = 0.064. 

As can be seen from figure 3, larger values of N s and £ both lead to a lower prediction 
for the primordial Li abundance. Taking for instance N s = 2, £ = 0.066 and uj\^ = 0.0225, 
one obtains [ 7 Li/H] p = 0.372 xlO" 10 (compared to [ 7 Li/H] p = 0.460 xlO" 10 ioi N s = £ = 0). 
While this is still 3.4 standard deviations larger than the measured value, the discrepancy is 
not quite as serious as the 4.8 standard deviations one finds in standard BBN. 
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3.6 Combining BBN with CMB+LSS data 

If one assumes that neither oj^ nor N e s changed between the BBN era and the time of 
photon decoupling, one could perform a combined analysis of the CMB+LSS and BBN data 
in analogy to the previous literature [44, 45]. 

With the stringent upper bound on the radiation density from the standard BBN setting, 
it is clear that the scenarios with one eV-mass sterile neutrino plus extra massless degrees of 
freedom considered in section 2 will fit the combined data rather badly: for a 1 eV (2 eV) 
sterile neutrino, the ACDM+AiV model yields AXg ff ^ 9 (19), relative to ACDM. This, in 
turn suggests the necessity for further modifications of the cosmological model: either by 
dropping the assumption of iV e g (or uj^) being constant, or, as seen above, by admitting a 
nonzero neutrino chemical potential. 

Aside of the massive sterile neutrino scenario, the combination of CMB+LSS and BBN 
data can of course also be used to constrain the commonly considered case of N e Q massless 
degrees of freedom. We find A^ cff = 3.90±g;|§ at 95% credibility, with the standard model 
expectation of N e g = 3.046 outside the 99.5%-credible region, similar to the limits reported 
in the recent work of Hou et al. [9] . 

4 Conclusions 

In this work we have investigated the effects of eV-mass sterile neutrinos, as suggested by 
global interpretations of neutrino oscillation data, on cosmology. Such sterile neutrinos can 
thermalise prior to neutrino decoupling, thus contributing to the relativistic energy density in 
the early universe. However, while the combination of CMB+LSS and BBN data does appear 
to prefer extra relativistic degrees of freedom at 99.5% credibility within the ACDM frame- 
work, fully thermalised massive sterile neutrinos in the 1-2 eV mass range necessarily violates 
the hot dark matter limit on the maximum neutrino mass. In terms of the goodness-of-fit, 
adding one massless sterile neutrino species improves the CMB+LSS fit by Ax^g = —3.16 
relative to standard ACDM, whereas endowing this sterile state with a mass of 1 eV worsens 
the fit by A^g ff = 4.20 relative to the same benchmark. Such a scenario is then excluded or 
strongly disfavoured. 

Nonetheless, while it appears difficult to accommodate eV-mass sterile neutrinos within 
the ACDM framework, extending the framework with modifications in the neutrino sector 
improves to some extent the consistency of sterile neutrinos with precision cosmological 
data. The simplest such modification is to admit even more additional (effectively massless) 
relativistic degrees of freedom, not necessarily fully thermalised. Such a scenario improves 
somewhat the bad effect of the 1 eV sterile neutrino mass at the expense of introducing an 
additional 1.5 massless species, but the xls °f the fit is still worse than standard ACDM by 
2.26 units. Allowing in addition for a dark energy equation of state w 7^ —1 further improves 
the fit: for a 1 eV sterile neutrino, a model with w = —1.11 and 1.23 additional massless 
species in fact fits the data marginally better than standard ACDM (A^g- = —0.78). 

Importantly, any model containing eV-mass sterile neutrinos will induce an upward shift 
in the cold dark matter density inferred from precision cosmological data. This shift can have 
important consequence for, e.g., the SUSY dark matter parameter space. 

However, as is well known, increasing the radiation content in the early universe can be 
problematic for BBN. We find that while standard BBN prefers roughly an extra relativistic 
degree of freedom, which we interpret here as a thermalised 1-2 eV sterile neutrino species, 
additional fully thermalised massless species are strongly disfavoured. Nevertheless, it is 
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possible to circumvent these BBN constraints with the introduction of a u e chemical potential, 
which itself could have been created by active-sterile oscillations in the early universe. Thus, 
the state of affairs can be summarised as follows. We need additional radiation to reduce the 
bad effect of hot dark matter on precision cosmology, and a small neutrino chemical potential 
to undo the bad effect of too much radiation on BBN. In principle, both of these ingredients 
can originate from the neutrino sector alone. 

In summary, it is not trivial to accommodate a strongly mixed eV-mass sterile neutrino 
in cosmology. Additional ingredients are required, such as additional radiation, a neutrino 
chemical potential, or a nontrivial w parameter. In all cases, significant changes in the 
inferred values of other a priori unrelated cosmological parameters are also incurred, e.g., an 
increase in the cold dark matter density. Thus, should the experimental indications for eV- 
mass sterile neutrinos become stronger, one must consider a fairly complex modification of the 
standard ACDM cosmology. On the observational side, the upcoming precision measurement 
of AiV m i by Planck [46, 47] remains one of the most promising windows to physics beyond 
the standard model. 
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